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By using a method, previously established to calculate electromagnetic fields, we compute the 
force of light upon a metallic particle. This procedure is based on both Maxwell's Stress Tensor and 
the Couple Dipole Method. With these tools, we study the force when the particle is over a flat 
dielectric surface. The multiple interaction of light between the particle and the surface is fully taken 
into account. The wave illuminating the particle is either evanescent or propagating depending an 
whether or not total internal reflection takes place. We analyze the behaviour of this force on either 
a small or a large particle in terms of the wavelength. A remarkable result obtained for evanescent 
field illumination, is that the force on a small silver particle can be either attractive or repulsive 
depending on the wavelength. This behaviour also varies as the particle becomes larger. 

PACS numbers: 78.70.-g, 41.20.-q, 42.50.Vk 



I. INTRODUCTION 



II. COMPUTATION OF THE OPTICAL FORCES 
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Since the first demonstration of particle manipulation 
by the action of optical forces, 1 ' 2 optical tweezers 3 and 
other configurations of light beams have been established 
to hold suspended particles like molecules , 4 or more re- 
cently, dielectric spheres. 5-7 Also, the possibilities of cre- 
ating microstructures by optical binding and resonance 
effects have been discussed 8-12 as well as the control of 
particles by evanescent waves. 13,14 Only a few works exist 
on the interpretation, prediction and control of the opti- 
cal force acting on a small particle on a plane surface. To 
our knowledge, the only theoretical works dealing with 
this subject are those of Refs. [15,16,17]. In Ref. [15] no 
multiple interaction of the light between the particle and 
the dielectric surface is considered. On the other hand, 
Ref. [16] deals with a 2-D situation. Only recently in 
Ref. [17] the full 3-D case with multiple scattering was 
addressed for dielectric particles. 

This work, extends the study of Ref. [17] to metallic 
particles and, as such, this is the first theoretical study 
of light action on a metallic particle. We shall there- 
fore present a rigorous procedure to evaluate the elec- 
tromagnetic force in three dimensions. Further, we shall 
analyze how this force depends on the wavelength, dis- 
tance between the particle and the surface, angle of in- 
cidence (whether the excitation is a plane propagating 
or an evanescent wave), and on the excitation of plas- 
mons on the sphere. We shall make use of the couple 
dipole method previously employed, and whose validity 
was analyzed in detail in Ref. [17]. 

In section II we introduce a brief outline on the method 
used to compute the optical force on a particle. We also 
write its expression from the dipole approximation for a 
metallic sphere in presence of a surface. Then, in Sec- 
tion III A we present the results and discussion obtained 
in the limit of a small sphere, and in Section IIIB we 
analyze the case of larger spheres compared to the wave- 
length. 



The Coupled Dipole Method (CDM) was introduced 
by Purcell and Pennypacker in 1973. 18 In this paper we 
use this procedure together with Maxwell's Stress Ten- 
sor (MST) 19 in order to compute the optical forces on a 
metallic object in the presence of a surface. Since we de- 
veloped this method in a previous paper, 17 we shall now 
outline only its main features. It should be remarked 
that all calculations next will be written in CGS units 
for an object in vacuum. 

The system under study is a sphere, represented by a 
cubic array of N polarizable subunits, above a dielectric 
flat surface. The field at each subunit can be written: 

N 

E(Ti,u) =E (ri,o;) + ^[S(r < ,r J -,a;) 

4=1 

+ T(r i ,r i ,w)]a i (w)E(r j ,w). (1) 

where Eo(i*j,cj) is the field at the position in the ab- 
sence of the scattering object, T and S are the field sus- 
ceptibilities associated to the free space 20 and the sur- 
face, 21,22 respectively. oii(u) is the polarizability of the 
i th subunit. Like in Ref. [17] we use the polarizability of 
the Clausius-Mossotti relation with the radiative reaction 
term given by Draine: 23 
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where ao holds the usual Clausius-Mossotti relation ao — 
a 3 (e — l)/(e + 2). 17,24 In a recent paper, 25 we have shown 
the importance to compute the optical forces taking into 
account the radiative reaction term in the equation for 
the polarizability of a sphere. For a metallic sphere, the 
polarizability is written as a = q:q(1 + (2/3)ifegQ!o)/D 
with D = 1 + (4/3)fc 3 3m(a ) + (4/9)/cg|a | 2 , where the 
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asterisk stands for the complex conjugate and 3m de- 
notes the imaginary part. 

The force 26 at each subunit is: 25 
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where k and I stand for the components along cither x, y 
or z, and 5fte denotes the real part. The object is a set 
of N small dipoles so that it is possible to compute the 
force on each one from Eq. (3). Hence, to obtain the total 
force on the particle, it suffices to sum the contributions 
F(rj) on each dipole. 



for s-polarization, with A = (1 — e)/(l + e) being the 
Fresnel coefficient of the surface. We have assumed a 
dielectric surface, hence A is real. 

From Fig. 1 is easy to see that k x is always real what- 
ever the angle 6, hence we can write Eqs. (4) and (6) for 
metallic particles as: 
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FIG. 1. Geometry of the configuration considered. Sphere 
of radius a on a dielectric flat surface (e = 2.25). The incident 
wave vector k is in the XZ plane. 

Being the object under study a small sphere located at 
ro = (0, 0, zq) (see Fig. I), we can employ the dipole ap- 
proximation, and hence use directly Eq. (3) with N=l. 
Within the static approximation for the field suscepti- 
bility associated to the surface (SAFSAS) (that is to 
say fco = 0), we have found an analytical expression for 
E(ro,w) that yields the force components: 17 
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for ,s-polarization. In Eq. (8) the factor in front of 
(feSm(ao) / (2D) + k x k^\ao\ 2 / (3D)) for the two polar- 
izations constitutes the field intensity at z . The first 
term within these square brackets corresponds to the ab- 
sorbing force whereas the second represents the scatter- 
ing force on the sphere. We see from Eqs. (8) and (9) 
that F x always has the sign of k x . Notice that it is not 
possible to write a general equation for the force along 
the Z-direction, as k z will be either real or imaginary, 
according to the angle of incidence. 



III. RESULTS AND DISCUSSION 



A. small particles 



We first address a small isolated silver particle with 
radius a =I0nm. In this case we can use the dipole ap- 
proximation, hence we consider Eqs. (4)-(7) with A = 0. 
Fig. 2a presents the polarizability modulus (|a |) of the 
sphere. The maximum of the curve corresponds to the 
plasmon resonance, i.e. when the dielectric constant is 
equal to -2 in Drudc's model. Notice that, in this model 
the dielectric constant is real, on using experimental val- 
ues, 27 the dielectric constant is complex and the reso- 
nance is not exactly at 3?e(e) = —2 but slightly shifted. 
In Fig. 2b we plot the real part of the polarizability 
(3?e(ao)), and Fig. 2c shows its imaginary part (SJm(ao)). 
Fig. 2d represents the force in free space computed from 
an exact Mie calculation (full line) and by the dipole 
approximation from Eqs. (4)- (7) with A = (dashed 
line) and Eq. (2) for a. In this case, the dipole approx- 
imation slightly departs from Mie's calculation between 
350nm and 375nm. We can compute the polarizability 
a from the first Mie coefficient a\ given by Dungey and 
Bohren (DB). 28 Therefore, the electric dipole polarizabil- 
ity is a — 3zai/(2fcg). 29 The symbol + in Fig. 2d corre- 
sponds to the DB polarizability and it is exactly coinci- 
dent with the Mie calculation. When the optical constant 
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of the metallic sphere is close to the plasmon resonance, 
the calculation from the Clausius-Mossotti relation with 
the radiative reaction term, departs from the exact cal- 
culation, even for a small radius. We shall next use the 
polarizability of DB. Analytical calculations will always 
be done with Eq. (2) to get simple expressions, and thus a 
better understanding of the physics involved. The curve 
of the force obtained by Mie's calculation has exactly the 
same shape as the imaginary part of the polarizability, 
this is due to the fact that for a small metallic sphere the 
absorbing force is larger than the scattering force. 
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FIG. 2. From top to bottom: the first three curves repre- 
sent the modulus, the real part, and the imaginary part of the 
polarizability of a silver sphere with radius a — lOnm versus 
the wavelength. The fourth curve is the force on this particle 
in free space. Plain line: Mie calculation, dashed line: po- 
larizability of Clausius-Mossotti relationm with the radiative 
reaction term, symbol +: DB polarizability. 
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FIG. 3. Force along the ^-direction on a silver sphere with 
a — lOnm versus distance z in the dipole approximation. The 
angle of incidence is 6 = 0°. With the static approximation 
(no symbol) and in an exact calculation (+) for S. Dashed line 
A =360nm, plain line A = 270nm, and dotted line A =500nm. 
The inset shows details of the zero force. 

Next, we consider the small sphere on a dielectric plane 
surface as shown by Fig. 1 . Illumination takes place from 
the dielectric side with 8 = 0°, hence in vacuum k z = ko 
and k x = 0. Fig. 3 represents the force in the Z-direction 
from Eq. (7) versus z for different wavelengths. Far from 
the surface the force tends to the Mie limit. Near the 
surface the force decreases, and, depending on the wave- 
length, it can become negative. For a better understand- 
ing of this force we write Eq. (7) as: 
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having made the approximation \a\ ~ \oco\, and D « 1 
since we have a small sphere compared to the wavelength 
(koa <C 1). The factor in front of the bracket corresponds 
to the intensity of the field at z - The first and the second 
terms within brackets represent the interaction between 
the dipole moment associated to the sphere and the in- 
cident field: the first term is the absorbing force whereas 
the second one corresponds to the scattering force, hence 
these forces are always positive. The third term is due 
to the interaction between the dipole and the field ra- 
diated by the dipole and reflected by the surface. We 
can consider this term as a gradient force exerted on the 
sphere due to itself via the surface. Hence, this force 
is always negative whatever the relative permittivity e. 
Since this term is proportional to 1/zq, it becomes more 
dominant as the sphere approaches the surface, hence the 
force decreases. To derive the point zq at which the force 
vanishes, if such a point exists, let us assume the scat- 
tering force smaller than the absorbing force, then from 
Eq. (10) the zero force is: 
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This equation always has a solution. We find z for the 
three wavelengths used to be: A = 270nm, z = 7.9nm, 
A = 360nm, z = 12.8nm, A = 500nm, z = 21.6nm. 
Now, zq (the location of the center of the sphere) must 
be larger than the radius a, or else, the sphere would be 
buried in the surface. Therefore the first of those values 
of Zo is not possible. Hence, the force is always positive. 
Thus, the distance between the sphere and the surface is 
z=2.8nm and 11.6nm for A =360nm and 500nm, respec- 
tively. These values are very close to those shown in the 
inset of Fig. 2. Notice that near the plasmon resonance 
for the sphere (A w 360nm) both |«o| and 9m («o) are 
maxima, hence the force is very large when the sphere 
is far from the surface and, due to the third term of 
Eq. (10), which depends of cto | 2 , the decay of this force 
is very fast. We have used the static approximation for 
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the susceptibility tensor associated to the surface (SAF- 
SAS) at large distance. We plot in Fig. 2 with crosses the 
force obtained from an exact calculation for S with the 
dipole approximation. As we see, these crosses coincide 
with those curves obtained with the SAFSAS whatever 
the distance z. Yet, we obtain for dielectric spheres the 
following: near the surface the SAFSAS is valid, whereas 
far from the surface the sphere does not feel its presence, 
and, thus, whether using the tensor susceptibility associ- 
ated to the surface in its exact form, or within the static 
approximation, has no influence. 



S 



o 

CD 

4 







J In 


















s — ■ 

I / — -1 







2 -4 



-6 

250 350 450 550 650 
wavelength (nm) 

FIG. 4. Force along the Z-direction on a silver sphere with 
a — lOnm versus the wavelenght A in the dipole approxi- 
mation. The angle of incidence is 6 — 50°. Plain line: 
p-polarization; dashed line: s-polarization. 
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FIG. 5. Force along the X-direction on a silver sphere 
with a = lOnm versus the wavelenght A in the dipole ap- 
proximation. The angle of incidence is 9 — 50°. Plain line: 
p-polarization; dashed line: s-polarization. 

We next consider the surface illuminated at angle of 
incidence 9 larger than the critical angle: 9 = 50° > 



41.8° = 9 C . Now the transmitted electromagnetic wave 
above the surface is evanescent. We plot in Fig. 4 from 
Eqs. (4)- (7) the force on the sphere for the two polar- 
izations in the Z-dircction versus the wavelength, and 
in the A-dircction in Fig. 5, when the sphere is located 
at zo=30nm. In Fig. 4 we see that the force in the in- 
direction is also cither positive or negative. In a previous 
work 17 , we have observed that the force on a small dielec- 
tric sphere is always attractive when the sphere is located 
in an evanescent wave. This is no longer the case for a 
metallic sphere. To understand this difference, and as 
the two polarizations have the same behaviour, we take 
the analytical solution for F z with k z = i-f (7 > 0) for 
s-polarization. Then Eq. (7) can be written: 
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As shown by Fig. 3, when the sphere is located at 
zq =30nm, the influence of the surface becomes negli- 
gible, thus we can use Eq. (13) with the hypothesis that 
z is large. Hence, F z w — \E 0y \ 2 -f$le(a )/2. This is 
the gradient force due to the incident field, and there- 
fore due to the interaction between the dipole associated 
to the sphere and the applied field. This force exactly 
follows the behaviour of 3?e(a ) (cf. Fig. 2). When 
5Re(«o) is negative, the dipole moment of the sphere 
oscillates in opposition to the applied field and so the 
sphere is attracted towards the weaker field. Notice that 
the same phenomenon is used to build an atomic mir- 
ror: for frequencies of oscillation higher than the atomic 
frequency of resonance, the induced dipole oscillates in 
phase opposition with respect to the field. The atom 
then undergoes a force directed towards the region of 
weaker field. 30 For p-polarization, the force can be writ- 
ten F z « -(l^oxl 2 + \E 0z \ 2 )-/^e(a )/2. As the modulus 
of the field becomes more predominant in p-polarization, 
the magnitude of the force becomes more important. We 
now search more carefully the change of sign in the force. 
Writing e = e' + is" for the relative permittivity, we get: 
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If the damping is weak, then the change of sign of F z 
happens both for e' « 1 and at the plasmon resonance 
for the sphere, i.e., e' w —2. Between these two values, 
the gradient force is positive. In fact, the limiting val- 
ues of the positive gradient force are always strictly in 
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the interval [—2, 1] due to damping. For example, the 
force vanishes at A = 352nm with e = —1.91 + 0.6i and 
A = 317nm with e = 0.66 + 0.95i. We notice that the 
change of sign happens steeply at the plasmon resonance 
since then the denominator of the real part of the po- 
larizability becomes very weak (see Eq. (14)), hence the 
zero force is surrounded by the two maxima of the force 
(one positive and the other negative). At A = 317nm 
the change of sign is smoother as in that case the de- 
nominator is far from zero. The third case, A = 259nm, 
lies between those two previous cases as the damping of 
the relative permittivity is important: e = —1.65 + 1.12«. 
We have also investigated the cases of gold and copper 
spheres, where a plasmon easily takes place, but we found 
not possible change of sign for F z for these two materials 
as the damping is then too important: if e" > 3/2, then 
Sfte(ao) is always positive whatever e' . However, if the 
particle is embedded in a liquid with a relative permit- 
tivity 2, then it is possible to get 5fte(a) < for gold. 
Notice that if 9 is close to 9 C , then 7 w and so we only 
have the third term of Eq. (13), then the force is always 
negative whatever the wavelength. In Fig. 5 we see, as 
previously, that the force in the A-direction has the sign 
of k x and, as the absorbing force is the most predomi- 
nant one, the curve has the same shape as the imaginary 
part of the polarizability (cf. Fig. 2). In that case, the 
maximum of the force F x is at the plasmon resonance 
(see Eq. (15)) namely at A = 354nm. 
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FIG. 6. Force along the Z-direction on a silver sphere 
with a — lOnm versus distance 2 in the dipole approxima- 
tion. The angle of incidence is 9 — 50°. Without symbol: 
p-polarization; with cross (+): s-polarization. Plain line: 
A =340nm, dashed line: A260=nm, dotted line: A =317. 5nm. 

Fig. 6 shows the force in the Z-direction versus z 
for both s (symbol +) and p-polarization (without any 
symbol) at three different wavelengths (A =265, 340, 
317. 5nm) for 9 = 50°. The behaviour is the same for 
both polarizations, only appearing a difference of magni- 
tude, this is due to the component of the field perpendic- 
ular to the surface in p-polarization. All curves manifest 



that near the surface the force is attractive, this is due, 
as seen before with a propagating wave, to the term of 
A/zq in Eqs. (10) and (13) which is always negative, 
irrespective of the kind of wave above the surface. At 
A =317. 5nm , (dotted line), we have 5Re(a) = which is 
why the force very quickly goes to zero when z grows. The 
two other cases correspond to 3?e(a) > (A =260nm) and 
Sfte(a) < (A =340nm) and far from the surface the force 
tends to zero. As the force is proportional to |E | 2 , we 
have F z cx e^ z \ 



B. large particles 

It is difficult to obtain convergence for the CDM when 
the relative permittivity of the medium to discretize is 
large. This imposes a very fine sampling. In this section, 
we use the range 250-355nm for the wavelength, the real 
part of the relative permittivity being small. Then, the 
computation is simplified. In addition, this is the most 
interesting case as then 5fte(a) = three times in this 
interval of wavelengths. 
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FIG. 7. Force along the Z-direction on a silver sphere ver- 
sus distance z with 9 = and a=100nm for the following 
wavelengths: A = 255nm (plain line), A = 300nm (dashed 
line), and A = 340nm (dotted line). Dots correspond to com- 
puted points. 

In Fig. 7 we plot the force in the Z-direction for an in- 
cident propagating wave (9 = 0°) versus the distance be- 
tween the sphere and the surface at three different wave- 
lengths: A =255nm, 300nm, 340nm. The calculations are 
done without any approximation. The curves have a sim- 
ilar magnitude and behaviour at the three wavelengths. 
The forces present oscillations due to the multiple reflec- 
tion of the radiative waves between the sphere and the 
surface, hence the period of these oscillations is A/2. The 
magnitude of these oscillations depends on the reflectiv- 
ity of the sphere, so the higher the Fresnel coefficient is, 
the longer these oscillations are. As expected, they are 
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less remarkable when the sphere goes far from the surface. 
We notice that the decay of the force when the metallic 
sphere gets close to the surface is not comparable to that 
on a dielectric sphere (see Ref.[17]). This is due to strong 
absorbing and scattering forces on the metallic sphere in 
comparison to the gradient force induced by the presence 
of the dielectric plane. 
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FIG. 8. Force along the Z-direction on a silver sphere with 
a=100nm versus distance z with 9 — 50° for the following 
wavelengths: A = 255nm (plain line), A = 300nm (dashed 
line), and A = 340nm (dotted line), symbol +: s-polarization; 
without symbol: p-polarization. 

In Fig. 8 we plot for 9 = 50° the normalized force 
in the Z-dircction, i.e., i^ z /|E | 2 , E being the field at 
zo in the absence of the sphere. We remark two im- 
portant facts at this angle of incidence. First, the de- 
cay of the force when the sphere is near the surface is 
more important in p-polarization. Notice that with the 
CDM it is not possible to numerically split the scatter- 
ing, absorbing, and gradient forces. Therefore, when the 
sphere is large we shall argument on the set of dipolcs 
forming it. In p-polarization, due to the ^-component 
of the incident field, the dipoles also have a component 
perpendicular to the surface, and this is larger than in 
s-polarization. Hence, in agreement with Fig. 6, due to 
this z-component, the attraction of the sphere towards 
the surface is larger for p-polarization. Second, all forces 
are positive when the sphere is far from the surface ex- 
cept for A =300nm in p-polarization. At this wavelength, 
for a small sphere, the force is negative for both s and 
p-polarization, hence we can assume an effect due to the 
size of the sphere. Just to see this effect, we present in 
Fig. 9 the force in the Z-direction versus the radius a, on 
a sphere located at zq — lOOnm for an angle of incidence 
9 = 50° but without taking into account the multiple 
interaction with the surface (i.e., S = 0). We take the 
previous wavelength of Fig. 8 (A = 255nm, 300nm, and 
340nm) more the wavelength at the plasmon resonance, 
A = 351. 5nm where 5?(ao) = 0. For small radius, we 



observe the same behaviour as in the previous section for 
an incident evanescent wave, namely at A = 255nm and 
340nm, the force in the Z-direction is positive, and it 
becomes negative for A = 300nm for both polarizations. 
These curves show a dependence proportional to the cube 
of the radius as 5fc(a ) oc a 3 . At the plasmon resonance, 
the force is slightly positive as there is no gradient force, 
but only weak absorbing and scattering forces. 
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FIG. 9. Force along the Z-direction on a silver sphere lo- 
cated at «o =100nm, with 8 = 50°, versus the radius a for 
A = 255nm (plain line), A = 300nm (dashed line), A = 340nm 
(dotted line), and A = 351. 5nm (thick line). Symbol +: 
s-polarization; without symbol: p-polarization. The inter- 
action between the sphere and the surface is not taken into 
account. 

But as shown by Fig. 9, when the radius grows, in s- 
polarization at A =300nm, the force sign changes and 
it becomes positive around 82nm as p-polarization keep 
the same behaviour. This confirms the fact that the pos- 
itive force obtained in Fig. 8 for A =300nm is only a size 
effect. For the cases A = 255nm, 340nm, the gradient 
force is positive in the Z-dircction, like the scattering 
and absorbing forces . Hence, the force, is always posi- 
tive whatever the radius. Nevertheless, for A = 300nm, 
these is a negative gradient force, the two other forces 
being positive. As previously said, it is not possible to 
know the relative contribution of the different forces, but 
in s-polarization the dipoles have a component mainly 
parallel to the plane. In that case, the dipole associated 
to each subunit radiates a field both propagating and 
evanescent along the Z-direction. Hence, we can assume 
that due to the radiative part, the absorbing and scatter- 
ing forces acting on the subunits inside the sphere, be- 
come important when the radius increases. Concerning 
p-polarization, due to the component of the field perpen- 
dicular to the plane, the dipoles have a component along 
Z and thus only radiate evanescent waves along this di- 
rection. In this case, we believe that the gradient force 
is more important on each subunit and counterbalances 
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the absorbing and scattering forces. At plasmon reso- 
nance A = 351. 5nm, when the radius grows, the absorb- 
ing and scattering forces become larger, but as no gra- 
dient force exists, the force is lower than those obtained 
at A = 255nm and 340nm. In fact, the curve shown at 
A = 340nm is the most important one, due to the onset 
of the plasmon resonance, thus |ao| and Sm(«o) being 
near their maximum, and 5Re(ao) close to its minimum. 
In this case, the gradient force is maximum and positive. 



IV. CONCLUSIONS 

We have presented a theoretical study of the optical 
forces acting upon a metallic particle on a dielectric plane 
surface cither illuminated at normal incidence or under 
total internal reflection. This study is done both with the 
coupled dipole method and Maxwell's stress tensor. We 
observe that when the incident wave is propagating, the 
difference between the force acting on a dielectric sphere 
and that on a metallic sphere stems from the absorbing 
force. Due to this contribution, the force upon a small 
silver sphere close to the dielectric surface can be posi- 
tive in spite of the gradient force. The opposite happens 
with a dielectric sphere. The main difference between 
the two cases (dielectric and metallic) arises however on 
illumination under total internal reflection. In that case, 
the effect on a small silver sphere is completely different 
to that observed on a dielectric sphere. Depending on 
the wavelength, the gradient force due to the incident 
field can be either repulsive or attractive. The change of 
sign happens both at the plasmon resonance and when 
e becomes close to one. In the interval between these 
two values the gradient force is positive. The explana- 
tion is very similar to that on the effect used to build 
an atomic mirror. At a wavelength where the gradient 
force on a small sphere is negative, we see that when 
the sphere radius grows, the force along the Z-direction 
stays negative for p-polarization, and it becomes positive 
for s-polarization due to the size effect. Nevertheless, at 
any arbitrary wavelength and angle of incidence, as the 
sphere approaches close to the surface, the attraction of 
the surface on the sphere increases: repulsive forces di- 
minish and even can change sign, eventually becoming 
attractive at certain wavelengths. Attractive forces, on 
the other hand, increase their magnitude. 
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